GL-TR-90-0349 


AD-A237  675 


Effects  of  Solar  Heating  by  Aerosols 

and  Trace  Gases  on  the  Temperature  Structure  Constant 


J.R.  Hummel 
E.P.  Shettle 


SPARTA,  Inc. 

24  Hartwell  Avenue 
Lexington,  MA  02173 


9  August  1990 


DTIC 

ELECTE 

JUL0  8199!. 


Scientific  Report  No,  6 


Approved  for  Public  Release;  Distribution  Unlimited 


GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 

HANSCOM  AIR  FORCE  BASE,  MASSACHUSETTS  01731-5000 


SI-03945 


This  technical  report  has  been  reviewed  and  approved  for  publication 


Contract  Manager  Electro-Optical  Measurement  Branch 


FOR  THE  COMMANDER 


R.  Earl  Good,  SES,  Director 
Optical  Environment  Division 


This  document  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) . 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center.  All  others  should  apply  to  the  National  Technical 
Information  Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  PL/IMA,  Hanscom  AFB,  MA  01731-5000.  This  will  assist  us  in  maintaining 
a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  L.e  .eturueu. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


1«.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBEfl(S) 

LTR90-009 


ea  NAME  OF  PERFORMING  ORGANIZATION 

SPARTA,  INC. 


6c.  ADDRESS  (City,  Suit,  and  Zip  Cod0) 

24  HARTAA'ELL  AVENUE 
LEXINGTON,  MA  02173 


Sa.  NAME  OF  FUNDING  /  SPONSORING  ORGANIZATION  8 


3c.  ADDRESS  (Ctty,  SUtt,  Mnd  Zip  Cod0) 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVElllARKINGS 


(It  0ppUe0bl0) 


(It  0ppllC0bl0) 


3.  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

GL-TR-90-0349 


7a.  NAME  OF  MONITORING  ORGANIZATION 

GEOPHYSICS  LABORATORY 


7b.  ADDRESS  (Clly,  St0l0,  0nd  Zip  Cod0) 

HANSCOM  AFB,  MA  01731-5000 


0.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

!  F19628-88-C-0038 


10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO. 

NO. 

ACCESSION  NO. 

621  OIF 

7670 

15 

AP 

11.  TITLE  (Includ0  S0cuniy  Cl000ilic0tianl 

EFFECTS  OF  SOLAR  HEATING  RV  AEROSOLS  AND  TRACE  GASES  ON  THE  TEMPERATURE  STRUCTURE 
CONSTANT 


12.  PERSONAL  AUTHOR(S) 

J  R  Hummel  and  E.  P.  Shettle* 


TIME  COVERED  14.  DATE  OF  REPORT  (Yaar,  Moom,  Day;  15.  PAGE  COUNT 

SCIENTIFIC  REPORT  #6  |  from  07/01/88  to  12/31/90  1990  AUGUST  9  24 

16.  SUPPLEMENTARY  NOTATION 

'Naval  Research  Laboratory,  Code  65221 
Washington  DC  20375 


1 8  SUBJECT  TERMS  {Conttnu0  cn  t1  n^cmtmmry  and  id^nufy  by  biock  number} 

■  '  10  I  GROUP  I  SUB  GROUP  | 

Optical  Turbulence,  Solar  Heating,  Ozone,  Aerosols,  Index  of  Refraction, 
Stnjcture  Function 


1^3.  ABSTRACT  {Cont$nu»  on  rovrsm  li  n0C9»ssry  snd  tdsntify  by  block  numbmrt 

j  ^■"hermosonde  measurements  made  in  the  stratosphere  indicate  an  increase  in  the  temperature  structure  constant, 
j  0  T  during  daylight  hours.  Simplified  calculations  indicate  that  temperature  perturbations  resulting  from  solar  heating  by 
jlhiii  layers  of  areosols  and  ozone  could  provide  the  mechanism  for  the  observed  day-night  differences  m  c  y 

Estimates  of  the  vertical  variations  in  ozone  and  aerosol  abundances  on  spatial  scales  between  a  few  tens  of  meters 
to  about  a  kilometer  have  been  used  with  mean  ozone  and  aerosol  profiles  to  make  detailed  calculations  of  the  solar 
heating  rates  by  ozone  and  aerosols.  The  calculations  iiave  utilized  the  temperature  dependent  ozone  absorption  cross 
sections  and  the  aerosol  properties  from  the  AFGL  model  atmospheres.  Perturbations  of  various  amounts  and  at 
d  fferent  attitudes  have  been  applied  to  the  mean  ozone  and  aerosol  profiles  to  study  the  impact  on  the  heating  rates  with 
rr.e  goal  of  trying  to  understand  the  stratospheric  thermosonde  measurements. 


uO  FORM  1473,  84  MAR 


83  APR  edition  may  bo  uGod  until  oxhaijotod 
All  other  editions  are  obsolete. 


GECURtTY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


Contents 


1.  INTRODUCTION  1 

1.1  Background  and  Purpose  of  the  Study  1 

1 .2  Organization  of  Report  1 

2.  BACKGROUND  AND  RELEVANT  OBSERVATIONS  2 

2.1  Background  2 

2.2  Relevant  Observations  3 

2.2.1  Ozone  Measurements  3 

2.2.2  Aerosol  Measurments  4 

2.3  Related  Studies  6 

3.  DETAILS  OF  CALCULATIONS  6 

4.  RESULTS  8 

4.1  Unperturbed  Conditions  8 

4.2  Perturbed  Conditions  8 

4.3  Determining  the  Solar  Heating  Rates  Required  to  Initiate  Turbulence  10 

4.4  Sensitivity  of  Results  to  Model  Calculations  1 1 

4.4.1  Impact  of  Local  Time  1 1 

4.4.2  Impact  of  Latitude  1 1 

4.4.3  Impact  of  Height  of  Perturbation  1 3 

4.4.4  Impact  of  Choice  of  Stratospheric  Aerosol  13 


111 


5.  SUMMARY  AND  CONCLUSIONS  14 

5.1  Summary  14 

5.2  Conclusions  14 

5.2.1  Role  of  Perturbations  in  Ozone  14 

5.2.2  Role  of  Perturbations  in  Stratospheric  Aerosols  14 

References  16 


IV 


1 .  An  Example  of  the  Daytime  Enhancement  in  as  Measured 
With  an  AFGL  Thermosonde 

2.  Balloon  Measurements  of  Ozone  Mixing  Ratio  Over 
(a.)  15  to  30  km  and  (b.)  24  -  26  km  From  March  1984 

3.  Dustsonde  Measurements  of  Aerosol  Concentration  as  a 
Function  of  Altitude 

4.  Calculated  Solar  Heating  Rates  for  the  Background,  Unperturbed 
Model  Atmosphere 

5.  Heating  Rate  as  a  Function  of  Altitude  Assuming  a  50  % 

Change  in  (a.)  the  Aerosol  and  (b.)  Ozone  Concentrations 

6.  Values  of  the  Critical  Temperature  Gradient  as  a  Function  of 
Wind  Shear  for  Three  Representative  Stratospheric  Temperatures 

7.  An  Example  of  Vertical  Wind  Shears  Measured  During 
the  AFGL  -  Penn  State  Measurement  Program 

8.  Heating  Rate  Gradient  as  a  Function  of  Local  Time  for  Two 
Different  Solar  Albedo  Conditions 


Acknowledgement 


The  authors  would  like  to  acknowledge  the  contributions  of  Mr.  James  H. 
Brown  and  Ms.  Gail  P.  Anderson  of  the  Optical  Environment  Division.  They 
provided  significant  suggestions  and  insight  to  the  research  and  their  coinributions 
and  time  are  greatly  appreciated. 


VI 


EFFECTS  OF  SOLAR  HEATING  BY  AEROSOLS  AND  TRACE 
GASES  ON  THE  TEMPERATURE  STRUCTURE  CONSTANT 

1.  INTRODUCTION 

1.1  Background  and  Purpose  of  the  Study 

The  propagation  of  radiation  through  the  atmosphere  can  be  affected  by  vari¬ 
ations  in  the  index  of  refraction.  These  variations  can  be  studied  by  measuring 
the  index  of  refraction  structure  constant,  C\.  The  index  of  refraction  structure 
constant  is  a  measure  of  the  variations  in  the  atmospheric  index  of  refraction  as  a 
function  of  altitude  and,  in  turn,  is  a  measure  of  the  amount  of  optical  turbulence 
present  in  the  atmosphere. 

Thermosonde  measurements  of  made  during  daylight  conditions  indicate  an 
enhancement  in  the  temperature  structure  constant  relative  to  measurements  made 
at  night.  Simplified  calculations  indicate  that  temperature  perturbations  resulting 
from  solar  heating  by  aerosols  and  ozone  could  provide  the  mechanism  for  the 
observed  day-night  differences  in  C7^.  The  purpose  of  this  study  is  to  perform 
detailed  calculations  of  the  solar  heating  rates  of  ozone  and  aerosols  on  spatial  and 
temporal  scales  appropriate  to  the  observed  variations  in  the  temperature  structure 
constant. 

1.2  Organization  of  Report 

Chapter  2  sets  the  background  for  the  phenomenon,  presents  some  relevant 
observations,  and  discusses  related  stucies  from  the  scientific  literanire.  Chapter  3 


] 


describes  the  methodology  used  to  perform  the  calculations  and  Chapter  4  discusses 
the  results  from  the  calculations.  Finally,  Chapter  5  provides  the  conclusions  and 
recommendations  for  additional  work. 

2.  BACKGROUND  AND  RELEVANT  OBSERVATIONS 
2.1  Background 

The  Air  Force  Geophysics  Laboratory  (AFGL)  has  undertaken  an  extensive 
measurement  program  of  over  the  years  utilizing  balloon-borne  thermosondes. 
Thermosondes  measure  the  temperature  difference  over  a  one  meter  horizontal 
path.^  The  mean  square  temperature  fluctuation  is  evaluated  and  over  the  one 
meter  distance  equals  the  temperature  structure  constant,  Cj.  The  temperature 
structure  constant  is  then  related  to  the  index  of  refraction  structure  constant  via 
the  relationship,  neglecting  the  small  contribution  from  H2O, 


where  2  is  the  altitude  in  km,  P  is  the  atmospheric  pressure  in  mb,  and  T  is  the 
atmospheric  temperature  in  K.  The  thermosonde  unit  is  carried  piggyback  fashion 
with  a  conventional  radiosonde  which  measures  the  pressure,  temperature,  relative 
humidity,  and  wind  speed  and  direction. 

The  AFGL  program  has  measured  at  different  latitudes  and  under  day  and 
night  conditions.  The  thermosonde  measurements  of  have  indicated  a  dramatic 
difference  between  the  averaged  day  and  night  conditions.  Between  about  15 
and  30  km,  mere  is  an  eniiancemenl  of  day  time  values  that  is  not  seen  in 

the  nighttime  values.  Figure  1  demonstrates  the  daytime  enhancement  in  by 
comparing  two  thermosonde  flights  taken  during  a  measurement  program  conducted 
at  Penn  State  in  1986.  The  two  thermosonde  missions  were  taken  on  the  same  day 
and  under  tlie  same  synoptic  weather  conditions. 

This  enhancement  is  believed  to  be  real  and  not  due  to  any  instrumental  effects. 
One  possible  explanation  for  the  daytime  enhancement  is  from  solar  heating  by 
atmospheric  gases,  such  as  ozone,  and/or  aerosols.  A  detailed  examination  of 
possible  instrumental  sources  of  error  has  not  indicated  an  instrumental  source  for 
the  observed  phenomenon.'^  GL  is  continuing  to  study  the  possibility  of  balloon 
thermal  wake  effects  on  the  turbulence  measurements.  A  series  of  dropsonde  tests 

Brown,  J.  H.,  Good,  R.  E.,  Bench,  P.  M.,  and  Faucher,  G.,  (1982)  “Sonde  Experiments  for  Com¬ 
parative  Measurements  of  Optical  Turbulence”,  AFGL-TR-82-(X)79,  ADA  118740,  24  February 
1982. 

Brown,  J.  H.,  Dewan,  E.,  Thomas,  P.,  and  Murphy,  E.,  (1989)  “Study  of  Possible  Solar  Headng 
Effects  on  Thermosonde  Probes  -  Error  Analysis”,  Geophysics  Laboratory,  Hanscom  AFB,  MA, 
GL-TR-89-0178,  10  July  1989,  ADA  218116. 


arc  planned  for  lale  1990  to  study  if  balloon  wake  effects  could  be  contributing 
to  the  enhancement.  The  fact  that  tliis  phenomenon  is  observed  independently  of 
location  supports  the  belief  of  a  solar  heating-driven  mechanism. 


Figure  1 .  An  Example  of  the  Daytime  Enhancement  in  as  Measured  With  an 
AFGL  Thermosonde 


'Ll  Relevant  Ooservations 

Two  pieces  of  observational  evidence  support  the  hypothesis  that  solar  heating  is 
the  driving  force.  One,  the  ozone  mixing  ratio  and  solar  heating  rates  are  increasing 
in  this  altitude  range.  Also,  stratospheric  aerosols  peak  in  this  altitude  range. 
Second,  measurements  have  shown  that  the  prohles  of  ozone  and  stratospheric 
aerosols  exhibit  structure  on  spatial  scales  small  enough  to  possibly  give  rise  to 
the  thermal  variation  that  ultimately  gives  rise  to  the  variation  in  die  temperature 
stnicture  constant. 

2.2.1  Ozone  Measurements 

The  vertical  profile  of  ozone  is  routinely  measured  using  balloon-borne  elec¬ 
trochemical  ozonesondes,  rockelsunur^.  m  rim  u'l.'avioiet  'bsonUi^m  oi  ^tometers. 


and  satellites '  'I’ho^e  techniques  indicate  that  staicture  exists  in  the  vertical 
ozone  protile 


Meas!;r'-'”e'''';  trou;  ;*  ‘vdion  ■  ;-nc  tnei.  rc>ohi:hM  i-vstem  develi'pcd  by  re¬ 
searchers  at  Harvard  Uno.ersity  deiinxnsirate  that  tiie  ozone  mixing  ratio  vanes  on 
a  spatial  scale  consistent  with  that  of  the  (\'~  data.  Figures  2(a.)  and  (b.  )  show 
a  protiie  of  ozone  mass  mixing  ratio  measured  by  the  Harvard  instrument.  'Hie 
data  were  obtained  over  Palestine,  Texas  at  the  end  of  .March  1984.  Figure  2  (a.) 
show's  the  prohle  over  the  15  to  30  km  region  and  Figure  2  (b.)  show  s  the  stmcture 
in  the  prohle  over  a  sebxted  2  km  region.  Figure  2  (a.)  also  shows  the  .M-'GL 
midiatitude  ozone  prohle  for  comparison  purposes.  '  The  data  shown  in  the  Figure 
exhibit  a  structure  in  the  ozone  mass  mixing  ratio  on  the  order  of  tens  of  meters. 

2.2.2  Aerosol  Measurement.s 


The  vertical  prohle  of  aerosols  in  the  stratosphere  is  routinely  measured  by 
lidar  systems,  balloon-borne  dustsondes.  and  satellites  (eg.')  These  measurement 
systems  also  demonstrate  tiiat  stmcture  exists  in  the  aerosol  protile.  In  the  case 
of  lidar  and  dustsonde  systems,  stmcture  on  the  order  of  lOO's  of  meters  can  be 
resolved.  Figure  3  shows  dustsonde  measurements  of  aerosol  concentration  in  the 
stratosphere  for  particles  with  radii  greater  than  0.15/im.’'  The  results  shown  in 
the  Figure  were  taken  prior  to  the  eruption  of  FI  Chichon  and  while  some  of  the 
stmcture  displayed  is  possibly  noise,  the  data  still  exhibit  the  hne  stmcture  that  is 
characteristic  of  aerosol  layers. 


■  Dutsch,  H.U.,  and  Ling.  ('h.  ( 1969)  Critical  Cotnpan.son  of  the  Detennination  of  Venical  Ozone 
Distnbution  by  the  I'mkelir  Method  and  by  the  Flecirochemical  Sonde,  Arm.  Geophy.s.,  25:211- 
214. 

’  McCormick,  .M.lb.  .Swissier,  T.J.,  i  lilsenratfi.  L..  Krueger.  and  Osborn.  .M.T.  t  1984)  Saiellite 
and  Correlative  .Measurements  of  Stratosphcnc  Ozone:  Compiuason  of  Measurements  .Made  by 
.SAGL,  ECC  Balloons.  Chemiluminescent,  and  Optical  Rocketsondes.  J.  Gcophys.  Res.,  89. 
5315-5.^20. 

'  Hilsenrath,  E.,  Attmannspachcr.  \V  ,  Bass,  A.,  Evans,  W..  Hagemeyer,  K.,  Barnes.  R..-\.,  Komhyr, 
W..  Mauersberger,  K..  .Mcntali,  .1.,  Rrofitt,  M.,  Robbins,  D..  Taylor,  S.,  Torres.  .A.,  and  WeinstcKk. 
E.,  0986)  Results  Fmm  the  Balloon  Ozone  Intercom  pan  son  Campaign  (BOIC),/  Gcophys.  Res.. 
91:13137-13152. 

Anderson,  G.P.,  Clough,  S.A..  Kneizys,  F.X.,  Cheiwynd,  J.H.,  and  Shettle,  E.P.,  (1986)  “.XFGL 
Atmospheric  Constilueni  ITotiics  d)-12()  km)".  .\l-CiE  rR-86-()l  10.  .ADA  175173,  15  May  1986 
'  Hummel,  J.R.,  .Sheltlc,  E.P  .  .md  l.ongim,  D.R.  ( 1988)  “A  .New  Background  Stratospheric  .Aerosol 
Model  for  Use  in  Atmospheric  Radiation  MihIcIs”.  Air  l-orce  Geophy''ics  Labc^ratory,  Hanscom 
AFB,  MA,  AFGL-TR-88  ()]66.  ADA  2101  10. 

^  Hofmann,  D.J.,  and  Rosen,  J  .\f,  (1983)  "Study  of  Large  .Aerosol  P.iniclcs'’,  Universitv  of 
Wyoming,  Laramie,  WY.  Report  «  AP-75. 
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(a.;  Ozone  Mixing  Ratio  (ppm) 


2  ^  Q  f  ^  t  I  i _  I  I  .  ‘  1  t  i  f  ‘  ■  I  '  ‘  ' _ 1  ;  ■  I  I.  A  i  i _ t  1  t  i  J _ 

4.5  5.0  5.5  6.0  6.5 


(b.)  Ozone  Mixing  Ratio  (ppm) 

Figure  2.  Balloon  Measurements  of  Ozone  Mixing  Ratio  Over  (a.)  15  to  30  km 
and  (b.)  24  -  26  km  From  March  1984.  (Data  Courtesy  of  Dr.  H.  Weinstock, 
Harv’ard  University.) 


Figure  3.  Dustsonde  Measurements  of  Aerosol  Concentration  as  a  Function  of 
Altitude 

2.3  Related  Studies 

Fiocco  et  have  studied  the  energy  exchange  that  occurs  with  aerosols 

in  the  atmosphere.  Their  studies  have  indicated  that  during  daylight  hours  aerosols 
can  achieve  temperatures  that  differ  from  the  ambient  atmosphere.  In  the  upper 
atmosphere,  the  solar  heated  aerosols  do  not  transfer  all  of  the  heat  energy  to  the 
ambient  atmosphere,  thereby  leading  to  the  large  aerosol-atmosphere  temperature 
differences.  For  altitudes  below  50  km,  the  aerosols  and  atmosphere  exchange  en¬ 
ergy  more  efficiently.  As  a  result,  the  difference  between  the  aerosol  temperatures 
and  the  ambient  atmosphere  is  less  than  a  degree  K  and  is  a  function  of  particle 
size.  The  larger  the  particle,  the  greater  the  temperature  difference. 


3.  DETAILS  OF  CALCULATIONS 

Calculations  were  performed  for  a  model  atmosphere  extending  from  the  surface 
to  35  km.  The  atmospheric  structure  was  based  on  the  AFGL  midlatitude  summer 
atmosphere.^  The  model  atmosphere  was  divided  into  139  layers.  The  layers 

Fiocco,  G.,  Grams,  G.,  and  Visconti,  V.,  (1975)  Equilibrium  Temperatures  of  Small  Particles  in 
the  Earth’s  Upper  Atmosphere  (50-110  km),  J.  Atmos.  Terr.  Phys.,  37:1327-1337 
Fiocco,  G.,  Grams,  G.,  and  Mugnai,  A.  (1976)  Energy  Exchange  and  Temperature  of  Aerosols 
in  the  Earth’s  Atmosphere  (0  -  60  km),  J.  Atmos.  Sci.,  33:2415-2424. 
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were  approximately  125  m  thick  between  15  -  30  km  and  about  lOOO  m  thick 
elsewhere.  Two  layers  were  combined  to  form  an  atmospheric  slab,  the  center  of 
which  is  where  the  heating  rates  were  calculated. 

Aerosols  were  introduced  at  all  altitudes.  The  type  and  amount  of  aerosol 
was  assumed  to  be  a  function  of  altitude.  'Fhe  aerosol  attenuation  propertie;'  were 
based  on  those  used  in  the  AFGL  model  atmospheres.^^ Between  0  -  2  km, 
rural  aerosols  based  on  a  relative  humidity  of  70  %  were  assumed.  In  the  free 
troposphere  (2-10  km),  tropospheric  aerosols  were  assumed.  Finally,  above  10 
km,  stratospheric  aerosols  were  assumed.  Spring/summer  conditions  were  assumed. 
The  type  of  aerosol,  however,  was  varied  in  order  to  study  the  sensitivity  of  the 
calculations  on  aerosol  attenuation  properties. 

A  total  of  106  spectral  intervals  between  0.135  -  0.755  fim  were  used  for  the 
ozone  and  molecular  oxygen  calculations.  'Fhe  .solar  lluxcs  were  based  on  tliose 
used  in  LOWTRAN  7''h  dlic  cross  sections  used  were  those  from  LOWTRAN  7 
and  FASCOD”'*’'  llie  ozone  cross  section  are  temperature  dependent'*'  although 
the  effect  of  the  temperature  dependence  is  negligible. 

The  aerosol  attenuation  was  calculated  over  15  wavelengths  between  0.2  - 
20  /an.  The  overlap  between  the  aerosols  and  gases  is  also  taken  into  account. 
For  the  results  discussed  in  this  paper,  the  calculations  where  performed  for  solar 
noon  conditions  at  equinox.  Also,  it  is  assumed  that  the  underlying  surface  has  a 
surface  albedo  of  30  %. 


^  Shettle,  E.  P..  and  Fenn,  R.W.  (1979)  '‘.Models  of  the  Aerosols  of  the  Lower  .^tmosphere  and 
the  Effects  of  Humidity  Vanations”,  AFGL-TR-79-0214,  ADA  08591 

Fenn.  R.W.,  Clough,  Gallery,  W.O.,  Good,  R.E..  Kneizys,  F.X.,  Mill,  J.D  ,  Rothman,  L.S., 
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4.  RESULTS 

4.1  Unperturbed  Conditions 

Figure  4  shows  the  calculated  solar  heating  rates  for  the  assumed  unperturbed 
atmospheric  conditions.  The  solid  curve  is  the  solar  heating  rate  from  ozone, 
molecular  oxygen  and  aerosols.  The  curve  marked  “Gases  only”  is  the  calculated 
solar  heating  rate  if  aerosols  were  not  present.  The  calculations  assumed  a  strato¬ 
spheric  aerosol  profile  appropriate  for  background  conditions  but  with  attenuation 
properties  representative  of  aged  volcanic  aerosols.  Aged  volcanic  aerosols  are 
considerably  more  absorbing  than  background  stratospheric  aerosols  and  these  re¬ 
sults  give  a  feel  for  the  relative  role  of  stratospheric  aerosols  in  determining  the 
stratospheric  solar  heating  rates.  These  calculations  demonstrate  the  apparent  in¬ 
sensitivity  of  the  total  solar  heating  rate  to  background  stratospheric  aerosols.  It 
should  be  stressed,  however,  that  these  calculations  are  representative  of  back¬ 
ground  stratospheric  conditions. 


MIDLATITUDE  SUMMER 


ffigure  4.  Calculated  Solar  Heating  Rates  for  the  Background,  Unperturbed  Model 
Atmosphere 

4.2  Perturbed  Conditions 

Figure  5  displays  heating  rates  for  the  model  atmosphere  with  perturbations  in 
the  (a.)  aerosol  and  (b.)  ozone  concentrations.  In  both  cases  it  is  assumed  that  the 
perturbation  occurs  over  a  20  m  layer  centered  at  25  km  and  that  the  peak-to-peak 
perturbation  is  50%  of  the  mean  concentrations. 
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Figure  5.  Heating  Rate  as  a  Function  of  Altitude  Assuming  a  50  %  Qiange  in  (a.) 
the  Aerosol  and  (b.)  Ozone  Concentrations 
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4.3  Determining  the  Solar  Heating  Rates  Required  to  Initiate  Turbulence 

Turbulence  is  initiated  when  the  local  Richardson  number,  Ri,  is  less  than  the 
critical  Richardson  number,  1/4.  The  local  Richardson  number  is  given  as  the 
ratio  of  the  bouyancy  to  shear  as 


g  is  the  acceleration  of  gravity,  T  is  the  temperature  in  K,  S  is  the  local  vertical 
wind  shear,  is  the  local  vertical  gradient  in  temperature,  and  7^  is  the  adiabatic 
lapse  rate,  0.0098  K/m.  For  instability  to  occur,  (1)  is  expressed  as 


Solving  (2)  for  the  temperature  gradient  required  for  instability  to  occur,  one  can 
solve  for  the  critical  temperature  gradient 


TS'^ 

‘^9 


For  representative  stratospheric  temperatures.  Figure  6  displays  the  values  of 
critical  lapse  rate  as  a  function  of  wind  shear.  Lapse  rates  to  the  right  of  the  curves 
represent  regions  of  instability  and  possible  turbulence  onset. 


WIND  SHEAR  (sec'') 

Figure  6.  Values  of  the  Critical  Temperature  Gradient  as  a  Function  of  Wind  Shear 
for  Three  Representative  Stratospheric  Temperatures 
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In  order  to  induce  turbulence,  the  layers  ot  ozone  and  aerosol  must  persist  long 
enough  to  produce  a  lapse  rate  over  the  layer  that  exceeds  the  critical  value.  Layers 
on  the  order  of  25  m  can  persist  for  as  long  as  12  to  43  minutes,  depending  upon 
the  degree  of  layer  spreading. 

If  an  aerosol  layer  20  m  thick  existed  witn  a  50  %  peak-to-peak  perturbation 
in  the  concentration,  the  resulting  difference  in  solar  heating  rate  (relative  to  the 
unperturbed  values)  would  be  1.5  x  lO"^  K  s~^  m~^  If  this  layer  persisted  for 
45  minutes,  the  ^  would  amount  to  4.0  x  10~^  K  m“^.  In  order  for  turbulence 
to  begin,  an  unrealistically  large  wind  shear  would  be  required.  This  result  indi¬ 
cates  that  it  is  unlikely  that  perturbations  in  the  background  stratospheric  aerosol 
abundances  are  responsible  for  the  enhancement  seen  in  C^. 

For  an  ozone  layer  20  m  thick  with  a  50  %  peak-to-peak  perturbation  in  the 
concentration,  the  resulting  difference  in  solar  heating  rate  is  -1.6  x  10“^  K  s“^ 
m“^.  If  this  layer  persisted  for  45  minutes,  the  ^  would  amount  to  -4.4  x  10~^ 
K  m-^ 

In  order  to  induce  turbulence  in  the  presence  of  the  above  perturbation,  wind 
shears  on  the  order  of  about  0.03  s~^  or  greater  would  be  required.  Stratospheric 
wind  shears  measured  during  the  Perm  State  program  routinely  exceeded  this  value, 
as  shown  in  Figure  7.  This  provides  support  to  the  hypothesis  that  solar  heating 
by  the  ozone  layers  could  be  inducing  turbulence. 

4.4  Sensitivity  of  Results  to  Model  Calculations 

4.4.1  Impact  of  Local  Time 

The  discussion  above  was  based  on  results  obtained  at  local  noon.  Figure  8 
presents  the  heating  rate  gradient,  a  function  of  local  time.  The  results 

are  given  for  two  surface  albedo  conditions,  an  ozone  perturbation  of  50%,  and  a 
midlatitude  atmosphere. 

The  heating  rate  gradient  decreases  as  the  local  time  increases  up  to  local  noon, 
after  which  it  increases.  The  variation  in  heating  rate  gradient  is  due  to  changes  in 
the  solar  zenith  angle  and  corresponding  changes  in  mass  path.  Such  variation  is 
consistent  with  changes  in  Cy,  with  solar  zenith  angle  as  discussed  in  Reference  2. 

4.4.2  Impact  of  Latitude 

The  discussion  above  was  based  on  results  obtained  for  a  latitude  of  45°  N. 
Calculations  were  performed  for  conditions  at  35°  N  and  yielded  small  increase 
in  the  heating  rate  gradient  of  16  %  relative  to  those  at  45°.  The  slight  increase 
is  due  to  the  larger  amounts  of  stratospheric  ozone  present  at  the  lower  latitude^" 
and  is  consistent  with  themiosonde  measurement  variations  with  latitude. 
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WIND  SHEAR  (sec"') 

Figure  7.  An  Example  of  Vertical  Wind  Shears  Measured  During  the  AFGL  - 
Penn  State  Measurement  Program.  'Fhe  data  shown  were  obtained  at  an  altitude 
resolution  of  about  250  m 


Local  Time 


Figure  8.  Heating  Rate  Gradient  as  a  Function  of  Local  Time  for  Two  Different 
Solar  Albedo  Conditions 
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4.4.3  Impact  of  the  Height  of  Perturbation 

The  height  where  the  perturbation  occurs  will  have  an  impact  on  the  results  of 
the  calculations.  This  is  due  to  the  fact  that  both  the  ozone  and  aerosol  profiles  are 
changing  in  the  stratosphere.  The  ozone  mixing  ratio  profile  is  increasing  in  the 
altitude  range  where  the  erihancements  are  observed.  Also,  the  stratospheric 
extinction  coefficient  profile  peaks  in  the  10  -  20  km  region,  the  peak  altitude 
depending  on  the  degree  of  volcanic  activity. 

If  the  ozone  perturbation  was  assumed  to  occur  at  20  km  rather  than  25 
km,  the  resultant  heating  rate  gradient  is  lowered  by  over  a  factor  of  three  to 
-4.8  X  10~'  Ks“^  m~^.  If  the  aerosol  perturbation  occurred  at  20  km,  the 
increase  in  heating  rate  gradient  is  negligible. 

4.4.4  Impact  of  Choice  of  Stratospheric  Aerosol 

Calculations  were  done  using  an  aerosol  profile  representative  of  fresh  and  aged 
aerosols  under  extreme  volcanic  conditions.  The  calculations  were  performed 
with  an  assumed  50%  perturbation  over  a  20  m  interval  centered  at  20  and  25  km. 

With  fresh  volcanic  aerosols  and  the  perturbation  centered  at  20  km  the  heating 
rate  gradient  is  about  -3  x  10"®  K  s~^  m"^.  This  is  5  -  6  times  larger  than  the 
heating  rate  gradient  due  to  a  perturbation  in  the  ozone  abundance  at  this  altitude. 
With  the  perturbation  centered  at  25  km  the  heating  rate  is  -3  x  10"'  K  s~^  m'  ^ 

When  aged  aerosols  are  present,  the  resultant  heating  rate  gradient  for  a 
50%  perturbation  in  the  aerosols  over  a  20  m  interval  centered  on  25  km  was 
-3.4  X  10~^  K  s"^  m-1.  This  is  over  twice  that  used  for  the  default  condi¬ 
tions  but  still  significantly  smaller  than  for  ozone  perturbations.  If  the  perturbation 
occurred  at  about  20  km,  on  the  other  hand,  the  altitude  where  the  enhanced  vol¬ 
canic  extinction  profile  is  maximum,  the  resultant  heating  rate  gradient  would  be 
-2.8  X  10"®  K  s"^  m~^.  This  value  is  on  the  order  of  that  found  for  ozone 
perturbations. 

What  these  results  indicate  is  that  large  amounts  of  volcanic  aerosols  could 
lead  to  heating  rate  gradients  large  enough  to  produce  lapse  rates  that  could  exceed 
the  critical  lapse  rates  required  to  produce  instabilities.  Again,  this  is  consistent 
with  the  larger  diurnal  variations  of  C“  as  measured  by  the  thermosonde  in  Hawaii 
where  the  variations  were  larger  and  where  daytime  values  began  departing  from 
nighttime  values  at  lower  altitudes. 
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5.  SUMMARY  AND  CONCLUSIONS 

5.1  Summary 

Thermosonde  measurements  of  C“  made  during  daylight  conditions  indicate  an 
enhancement  in  the  temperature  structure  constant  relative  to  measurements  made 
at  night.  Simplified  calculations  indicate  that  temperature  perturbations  resulting 
from  solar  heating  by  aerosols  and  ozone  could  provide  the  mechanism  for  the 
observed  day-night  differences  in  However,  no  estimates  have  been  made 
for  C“  resulting  from  this  mechanism,  rather,  the  possibility  is  established  for  the 
generation  of  turbulence. 

'IXvo  pieces  of  observational  evidence  support  the  hypothesis  that  soiar  heating 
is  the  driving  force.  One,  the  ozone  mixing  ratio  and  solar  heating  rates  due  to 
ozone  are  increasing  in  this  altitude  range.  Also,  stratospheric  aerosols  peak  in  this 
altitude  range.  Second,  measurements  have  shown  that  the  profiles  of  ozone  and 
stratospheric  aerosols  exhibit  structure  on  spatial  scales  small  enough  to  possibly 
give  rise  to  the  thermal  variation  that  ultimately  gives  rise  to  the  variation  in  the 
temperature  structure  constant. 

5.2  Conclusions 

5.2.1  Role  of  Perturbations  in  Ozone 

The  results  for  assessing  the  role  of  perturbations  in  the  ozone  solar  heating 
are  encouraging.  Calculations  of  lapse  rates  that  would  result  from  solar  heating 
effects  alone  indicate  that  layers  on  the  order  of  20  m  with  perturbations  in  the 
ozone  profile  as  high  as  50  %  could  give  rise  to  turbulence  onset  in  the  stratosphere. 

5.2.2  Role  of  Perturbations  in  Stratospheric  Aerosols 

A  perturbation  of  50%  in  the  background  stratospheric  aerosol  concentration 
over  a  20  m  interval  centered  at  25  km  produces  a  change  in  solar  heating  rate  on 
the  order  of  -10~®  K  s~^  m''^,  two  orders  of  magnitude  smaller  than  that  from 
a  corresponding  50%  perturbation  in  ozone.  When  the  perturbation  occurs  at  20 
km,  the  heating  rate  gradient  is  still  on  the  order  of  -10“®  K  s"^  m“^  These 
results  indicate  that  it  is  unlikely  that  perturbations  in  the  background  stratospheric 
aerosols  could  be  responsible  for  the  observed  variations  in  the  temperature  struc¬ 
ture  constant. 

When  an  extreme  volcanic  aerosol  profile  is  assumed,  the  heating  rate  gradient 
is  on  the  order  of  -10“'  to  -10“®K  s“^  m“^,  depending  on  the  type  of  volcanic 
aerosols  present,  for  a  perturbation  occurring  at  25  km  but  is  on  the  order  of  - 
10“^  K  s“^  m“^  which  is  on  the  order  of  or  larger  than  that  from  perturbations 
in  the  ozone  abundance.  These  results  indicate  that  in  the  presence  of  extreme 
volcanic  conditions,  perturbations  in  the  aerosol  abundances  could  give  rise  to  the 


enhancement  in  C“. 

The  results  for  assessing  the  role  of  perturbations  in  the  ozone  solar  heating 
are  encouraging.  Calculations  of  lapse  rates  that  would  result  from  solar  heating 
effects  alone  indicate  that  layers  on  the  order  of  20  m  with  ozone  perturbations  as 
high  as  50  %  could  give  rise  to  turbulence  onset  in  the  stratosphere.  Thermosonde 
measurements  taken  in  Hawaii  at  a  time  when  the  Kilauea  volcano  was  active  are 
consistent  with  larger  diurnal  variations  beginning  at  lower  altitudes. 
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